The Isospin Dependence Of The Nuclear Equation Of State Near The Critical Point 
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We discuss experimental evidence for a nuclear phase transition driven by the different concentra- 
tion of neutrons to protons. Different ratios of the neutron to proton concentrations lead to different 
critical points for the phase transition. This is analogous to the phase transitions occurring in 4 He- 
3 He liquid mixtures. We present experimental results which reveal the N/A (or Z/A) dependence of 
the phase transition and discuss possible implications of these observations in terms of the Landau 
Free Energy description of critical phenomena. 
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Nuclei are quantum Fermi systems which exhibit ma- 
ny interesting features which depend on temperature 
and density. At zero temperature and ground state 
density, nuclei are charged quantum drops, i.e. they 
have a Fermi motion |l| due to their quantum nature, 
and nucleons interact through a short range attractive 
force and the long range Coulomb repulsions among the 
constituent protons. In the absence of the Coulomb 
force, the Nuclear Hamiltonian is perfectly symmetric 
for exchange of protons and neutrons apart from a small 
but NOT insignificant difference between the proton and 
neutron masses. This symmetry is revealed by similar 
energy levels in mirror nuclei, i.e. nuclei with the same 
mass number A but opposite numbers of neutrons, N, 
and protons, Z. Of course this feature is observed for 
relatively small systems because the Coulomb energy is 
small Analogous to the properties of mirror nuclei, 
we could expect that if we study nuclei at finite tem- 
peratures, T, and low densities, p, then, if the Coulomb 
force is not important, the invariance under exchange 
of protons to neutrons might lead to important and 
interesting consequences. In fact, since the fundamental 
Hamiltonian of nuclei is invariant under exchange of N 
with Z (apart from Coulomb effects), we could expect 
that such an invariance should be manifested only at 
high T (disordered state), while there is a spontaneous 
symmetry breaking at lower T (ordered state). That 
means that in symmetric nuclear matter at high T, 
the state with fragments having N — Z defines the 
minimum of the free energy, i.e. symmetric fragments 
such as deuterons and alphas would be favored at low 
density [H, H|. On the other hand, there could be a 
symmetry breaking favoring N ^= Z at lower T. In this 



case fragments near a (first order) phase transition might 
prefer either a neutron or proton rich configuration. 
There might even be a more interesting situation, 
suggested by the present data, the existence of a line 
of first order phase transitions [2| which terminates 
in a tri-critical point. For such a line the free energy 
has three equal minima: one with N=Z and the other 
two for N =/= Z. Thus a phase transition is driven by 
the difference in isospin concentration of the fragments 
to = (N — Z)/A. In this paper we will discuss data 
which clearly demonstrate that m is an order parameter 
of the phase transition. Its conjugate field Q which 
we indicate with H, is due to the chemical potential 
difference between protons and neutrons of the emitting 
source at the density and temperature reached during 
a collision between heavy ions 0, @ • We also note that 
the phase transition has a strong resemblance to that 
observed in superfluid mixtures of liquid 4 He- 3 He near 
the A point. In both systems, changing the concentration 
of one of the components of the mixture, changes the 
characteristics of the Equation Of State (EOS) 

In recent times a large body of experimental evidence 
has been interpreted as demonstrating the occurrence 
of a phase transition in finite nuclei at temperatures 
(T) of the order of 6 MeV and at densities, p, less 
than half of the normal ground state nuclear den- 
sity @. Even though strong signals for a first and a 
second-order phase transition have been found 0, 0], 
there remain a number of open questions regarding the 
Equation of State of nuclear matter near the critical 
point. In particular the roles of Coulomb, symmetry, 
pairing and shell effects have yet to be clearly delineated. 
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Theoretical modeling indicates that a nucleus excited 
in a collision expands nearly adiabatically until it is 
close to the instability region thus the expansion is 
isentropic [1] . At the last stage of the expansion the role 
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of the Coulomb force becomes very important. In fact, 
without the Coulomb force, the system would require 
a much larger initial compression and/or temperature 
in order to enter the instability region and fragment. 
The Coulomb force acts as an external piston, giving 
the system an 'extra push' to finally fragment. These 
features are clearly seen in Classical Molecular Dynam- 
ics (CMD) simulations of expanding drops with and 
without a Coulomb field 0, [l(| • The expansion with the 
Coulomb force included is very slow in the later stage 
and nearly isothermal. 

Even though at high T and small p the nucleus behaves 
as a classical fluid, the analogy to classical systems should 
not be overemphasized as, in the (T,p) region of interest, 
the nucleus is still a strongly interacting quantum sys- 
tem. In particular the ratio of T to the Fermi energy 
at the (presumed) critical point is still smaller than 1 
which suggests that the EOS of a nuclear system is quite 
different from the classical one. To date this expected 
difference has not been well explored @, [ill - fit^ . 

The paper is organized as follows: in the next sec- 
tion we discuss the experimental setup in detail. This is 
followed by a description of the data analysis and a dis- 
cussion in terms of the Landau 0(m 6 ) free energy. We, 
then derive some critical exponents and the EOS corre- 
sponding to possible scenarios suggested by our data in 
terms of the Fisher model of fragmentation. Finally we 
draw some conclusions and suggest possible future work. 



I. EXPERIMENTAL DETAILS 

The experiment was performed at the K-500 super- 
conducting cyclotron facility at Texas AhM Univer- 
sity. 64 ' 70 Zn and 64 Ni beams were incident on 58 - 64 Ni, 
ii2,i24 Snj i97 Au and 232 Th targets at 40 A MeV. Interme- 
diate mass fragments (IMF) were detected by a detector 
telescope placed at 20°. The telescope consisted of four 
Si detectors. Each Si detector had 5cm x 5cm area. The 
thicknesses were 129, 300, 1000 and 1000 fim. All Si de- 
tectors were segmented into four sections and each quad- 
rant had a 5° acceptance in polar and azimuthal angles. 
The fragments were detected at average angles of 17.5° 
± 2.5° and 22.5° ± 2.5°. Typically 6-8 isotopes were 
clearly identified for a given Z up to Z=18 with an en- 
ergy threshold of 4-10 A MeV, using the AE-E technique 
for any two consecutive detectors. The AE-E spectrum 
was linearized by an empirical code based on a range- 
energy table. In the code, isotopes are identified by a 
parameter Z^eai- For the isotopes with A=2Z, Zfieai = 
Z is assigned and other isotopes are identified by inter- 
polating between them. The energy spectrum of each 
isotope was extracted by gating on lines corresponding 
to the individual identified isotopes. In order to com- 
pensate for the imperfectness of the linearization, actual 
gates for isotopes were made on the 2D plot of Zn ea i 
versus energy. The multiplicity of each isotope was eval- 



uated from the extracted energy spectra, using a mov- 
ing source fit at the two given angles. Since the energy 
spectra of some isotopes have very low statistics, the fol- 
lowing procedure was adopted for the fits. Using a single 
source with a smeared source velocity around half of the 
beam velocity, the fit parameters were first determined 
from the energy spectrum summed over all isotopes for a 
given Z, assuming A=2Z. Then assuming that the shape 
of the velocity spectrum is the same for all isotopes for 
a given Z. All parameters except the normalizing multi- 
plicity parameter were assumed to be the same as for the 
summed spectrum. The multiplicity for a given isotope 
was then derived by normalizing the standard spectrum 
to the observed spectrum for that isotope. 

In order to evaluate the back ground contribution to 
the extracted multiplicity a two Gaussian fit to each iso- 
tope peak was used with a linear background. The second 
Gaussian (about 10% of the height of the first one) was 
added to reproduce the valleys between isotopes. This 
component was attributed to the reactions of the isotope 
in the Si detector. The centroid of the main Gaussian 
was set to the value calculated from the range-energy ta- 
ble within a small margin. The final multiplicity of an 
isotope with Z > 2 was obtained by correction of the 
multiplicity evaluated form the moving source fit for the 
ratio between the sum of the two Gaussian yields and the 
linear background. 

The yields of light charged particles (Z < 2) in coin- 
cidence with IMFs were also measured using 16 single 
crystal CsI(Tl) detectors of 3cm thickness set around the 
target. The light output from each detector was read 
by a photo multiplier tube. The pulse shape discrim- 
ination method was used to identify p, d, t, h and a 
particles. The energy calibration for these particles were 
performed using Si detectors (50 -300 /im) in front of the 
Csl detectors in separate runs. The yield of each iso- 
tope was evaluated, using a moving source fit. Three 
sources (projectile-like(PLF), nucleon-nucleon-like(NN) 
and target-like (TLF)) were used. The NN-like sources 
have source velocities of about a half of the beam ve- 
locity. The parameters were searched globally for all 16 
angles. Detailed procedures of the data analysis are also 
given in refs. [l7l [T^ |. 

Special care has been taken for 8 He identification. All 
He isotopes are identified in the Si telescope, using the 
AE-E technique, in a narrow energy range. When a pro- 
ton and an a hit the same quadrant and when both of 
them stop in the E detector, their AE-E points overlap 
with those of 8 He. Since the multiplicities of protons 
and alphas are about three orders of magnitude larger 
than that of 8 He, the contribution of accidental events 
becomes significant, especially for the reaction systems 
with lower numbers of neutrons, in which 8 He production 
is suppressed. Since Z=l AE-E spectra are not available 
in this experiment, AE-E for Z < 3 were measured in a 
separate run. Using the light charged particle multiplic- 
ity extracted from the 16 Csl detectors, the accidental 
events were simulated for each reaction for the observed 
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a yield in the AE-E spectra in this experiment, the solid 
angle of the quadrant and the multiplicity of Z—l parti- 
cles. In order to minimize the accidental events, the runs 
with a low beam intensity were selected in each reaction. 
Typical linearized Zft ea i spectra with these accidentals 
are shown in Fig.{TJ for 70 Zn + 232 Th (N/Z - 1.5) and 
64 Ni + 112 Sn ( N/Z - 1.25). As one can see, the Z rea i 
values for the accidental events of proton and a pileup 
is nearly identical to that of 8 He, while 6 He is clearly 
identified. The contributions from d + a and t 4- a are 
also reasonably consistent with the observed background 
yields. A significant excess of 8 He yield beyond the ac- 
cidentals is only observed for the reaction systems with 
the 124 Sn, 197 Au and 232 Th targets. After the correction 
of the accidental contributions, the multiplicities of 6 He 
and 8 He were calculated using the source fit parameters 
obtained for Li isotopes. 
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FIG. 1: Typical Zn ea i spectra for He and Li isotopes. Ac- 
cidental events are generated only for p, d, t +a and shown 
separately by shaded histograms as indicated. 



II. DATA ANALYSIS 

The key factor in our analysis is the value I = N — Z 
of the detected fragments. A plot of the yield versus 
mass number when 1=0 displays a power law behavior 
with yields decreasing as A~ T 0, This is shown in 
Fig.© for the 64 Ni+ 124 Sn case at 40 MeV/nucleon. In 
the figure we have made separate fits for odd-odd (open 
symbols) or even-even(filled symbols) nuclei. As seen, 
different exponents r appear which suggests that pairing 
is playing a role in the dynamics [l|, leading to higher 
yields for even-even nuclei. 

The observation of the power law behavior suggests 
that the mass distributions may be discussed in terms of 
a modified Fisher model @, [3 : 

Y = y A- r e~^ A , (1) 

where yo is a normalization constant, r = 2.3 is a crit- 
ical exponent Q, /3 is the inverse temperature and 



FIG. 2: Mass distribution for the 64 Ni+ 124 Sn system at 40 
MeV/nucleon for 1=0. The lines are power law fits with ex- 
ponents 2.3L0.02 (odd-odd nuclei, dashed line), and 3.4L0.06 
(even-even nuclei, full line) respectively. 



A/i = F(I/A) is the free energy per particle, F, near 
the critical point. Recall that in general, the free energy 
is a function of the mass A (volume), A 2 / 3 (surface) and 
the chemical composition m of the fragments and pos- 
sibly pairing. The region we are studying in this paper 
seems near the critical point for a liquid gas phase tran- 
sition (volume and surface equal to zero) but modified 
by m = Because of this modification we can observe 
different features of the transition such as a first order 
phase transition driven by m, the order parameter. 

We begin our analysis by noting that the Fisher free 
energy is usually written in terms of the volume and the 
surface of a drop undergoing a (second order) phase tran- 
sition 21] . Our data indicate that those terms are not im- 
portant in the present case [3] as we will show more in de- 
tail below. If they are negligible this suggests that we are 
near the critical point for a liquid gas phase transition. 
Because we have two different interacting fluids, neutrons 
and protons, the transition becomes more complex and 
more interesting than in a single component liquid. Ex- 
periments at different energies might display a free energy 
which depends on all these factors. If we accept that 
F is dominated by the symmetry energy we can make 
the approximation that F(I/A) — E sym — 25(I/A) 2 
MeV/A, i.e. the symmetry energy of a nucleus in its 
ground state [l|. We will use this relationship in order 
to infer an approximate value of the temperature of the 
system. However, we stress that in actuality, F(I/A) is 
a function of density, temperature and all other relevant 
quantities near the critical point. According to the Fisher 
equation given above, we can compare all systems on the 
same basis by normalizing the yields and factoring out 
the power law term. For this purpose we have chosen to 
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normalize the yield data for each system to the 12 C yield 
(I = 0) in that system, i.e. we define a ratio: 



R = 



Y{ 12 C)\2 T ' 



(2) 



The normalized ratios for the system 64 Ni + 64 Ni at 
40 MeV/nucleon are plotted as a function of the (ground 
state) symmetry energy in Fig.©, bottom panel. The 
data display an exponential decrease with increasing 
symmetry energy, except for the isotopes for which 1 = 0. 
The yields of these 1 = isotopes are of course not sensi- 
tive to the symmetry energy but rather to the Coulomb 
and pairing energies and possibly to shell effects. A fit 
to the exponentially decreasing portion of the data us- 
ing the ground state symmetry energy gives an 'appar- 
ent temperature' T of 6.0 MeV. This value of T would 
be the real one if only the symmetry energy is impor- 
tant, if entropy can be neglected, if a sym = 25 MeV (the 
g.s. symmetry energy coefficient value) and if secondary 
decay effects are negligible. In general we expect that 
the symmetry energy coefficient is density and tempera- 
ture dependent. Further, secondary decay processes may 
modify the primary fragment distributions (l7l [l8j . we 
will discuss these questions in the framework of the Lan- 
dau free energy approach below. We stress that the ap- 
pearance of two branches in FigJ3] (bottom), indicates 
that the total free energy must contain an odd power 
term in (I/A) at variance with the common expression 
for the ground state symmetry energy. For reference in 
the top part of Fig.© we have plotted the ratio versus 
the total ground state binding energy of the fragments. 
No clear correlations are observed which might suggest 
that the symmetry energy dominates the process. 
It is surprising that such a scaling appears as a function 
of the symmetry energy only. In fact we might won- 
der about the role of the Coulomb energy if we accept 
that surface and volume terms give negligible contribu- 
tion. In figure U we have plotted the same normalized 
ratios as a function of the quantity aE cou i + /3E sym , a 
and (5 are arbitrary parameters given in the figure and 
E cou i = 0.7Z(Z—1)A~ 1 ^ 3 is the Coulomb contribution to 
the ground state energy of the nucleus. We see from the 
figure that by decreasing the relative contribution of the 
Coulomb energy compared to the symmetry energy the 
scaling appears. This implies that the Coulomb energy 
is much less important than the symmetry energy near 
the critical point, which suggests that either the density 
dependence of those two terms is different or that, at the 
time of formation the fragments are strongly deformed, 
reducing the Coulomb effect. Such deformations have 
been seen in CMD calculations of fragmentation Q . 
To further explore the role of the relative nucleon con- 
centrations we plot in Fig.© the quantity ^ = — ln W 
versus m = {I /A), the difference in neutron and proton 
concentration of the fragment. As expected the normal- 
ized yield ratios depend strongly on m. 

Pursuing the question of phase transition we can per- 
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FIG. 3: Ratio versus fragments ground state binding energy 
(top panel) and symmetry energy (bottom panel) for the 64 Ni 
+ 64 Ni case at 40 MeV/nucleon. a s ym=25 MeV is used. The 
J < and I > (I = 0)isotopes are indicated by the open 
and full circles respectively (full squares). The dashed lines 
(bottom panel) are fits using a ground state symmetry energy, 
Eq[T] and a 'temperature' of 6 MeV. Notice that the given 
experimental 8 He yield is the upper limit. 



form a fit to these data within the generalized Landau 
free energy description [2j. In this approach the ratio of 
the free energy to the temperature is written in terms of 
an expansion: 



F 1 2 
— = —am 
T 2 
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bm 4 
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H 



— cm — TO—. 
6 T 



(3) 



where m is an order parameter, H is its conjugate vari- 
able and a — c are fitting parameters 0- We observe that 
the Free energy is even in the exchange of m — > —to, 
reflecting the invariance of the nuclear forces when ex- 
changing N and Z. This symmetry is violated by the 
conjugate field H which arises when the source is asym- 
metric in chemical composition. We stress that m and H 
are related to each other through the relation to = — 1§ . 

The use of the Landau approach is for guidance oifly. 
While the approximation to 0(m 4 ) does not work [I], 
the 0(m 6 ) case is in good agreement with the data. This 
is not surprising since, if fluctuations are important, a 
higher order approximation to the free energy is better, 
i.e. gives critical exponents closer to those seen in the 
data and satisfies the Ginzburg criterion [2]. A free fit 
using EqJ3] is displayed in FigJS] (full line). Notice the 
change of curvature near m = 0.3, which incidentally is 
close to m cn of the compound nucleus. For comparison 
in the same figure we have displayed the 0(m 2 ) case, i.e., 
F/T = a(m - m s ) 2 (b = c = 0) [26] As seen in the plot 
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FIG. 5: Free energy versus m for the case 64 Ni+ 232 Th. The 
full line is a free fit based on Landau 0(m 6 ) free energy. The 
dashed-dotted-dotted-dotted line is obtained imposing in the 
fit b = — ^/l6/3ac and it is located on a line of first order 
phase transitions. The short dashed line corresponds to b = 
— V4ac, i.e. superheating. The 0(m 2 ) case, F/T = a(m — 
m s ) 2 ,i.e. b = c = 0, m s — 0.1, is given by the long dashed 
line. 



ing values of m [H : 



last assumption also produces a reasonable fit, although 
it does not reproduce shoulders near m ~ ±0.3. As we 
will discuss in more detail below the appearance of two 
minima for m^0 (when H/T = 0) might be a signature 
for the existence of a first order phase transition occurring 
in these reactions. 

In general the coefficients entering the Landau free en- 
ergy Eq.©, depend on temperature, pressure or den- 
sity of the source. Usually one assumes c > 0, a = 
a<j(p)(T — To) and b = b(T, p), where To is some 'critical' 
temperature discussed below. The precise determination 
of these parameters determines the nuclear equation of 
state (NEOS) near the critical point. The data we have 
do not allow such a complete constraining of the NEOS 
but do suggest some interesting possible scenarios which 
we discuss below. 

We begin by noting that the conjugate variable H 
which appears in equation (J3|) is determined by the chem- 
ical composition of the source. Since, in general, the 
source has N ^ Z, the extreme of F/T are displaced 
from the values obtained when H=0. In fact if we take 
the first derivative of the free energy we get: 



(£)' 



bm + cm 



H 
T' 



(4) 



mo = 0; m ± = 



-b ± y/b 2 - 4ac 
2c ' 



(5) 



If we now assume H ^ but small, we can expand the 
solutions above as m = mo± + r\ with r\ small. Equating 
the first derivative to zero, Eq.Q, and neglecting terms 
0(tj 2 ) we get: 



V = 



H/T 



36mQ ± 



5crno ± 



(6) 



When H/T = the first derivative is zero for the follow- 



The shift of the minimum from uiq = should be given 
by the equation above and should be proportional to m of 
the emitting source. We can easily check this feature in 
our data. In Fig.© we plot the values of H/T obtained 
from the fits to our data for all systems using Eq. (j3]) ver- 
sus m cn = (I/A) cn . 

The linear fit in FigEO is given by H/T = 0.47 + 
1.6(I/A) cn which agrees with the linear dependence of 
Eq.©. However, for this fit H/T ^ for I cn = which 
could indicate the favoring of N > Z fragments by the 
Coulomb field. Another possibility is that (I / A) source oc 
(I/A)cn which then gives H = when I source = 0. 
Finally we should consider that together with H also 
the temperature may also be changing some since the 
collisions are between different target-projectile combi- 
nations at the same beam energy. If the temperature 
were the same then the coefficients of the free energy, 
EqJ21 should be independent of the source size, only H/T 
should change. In Fig. ([7]) we plot the parameters a, b 
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FIG. 6: H/T versus (I/A) of the compound nucleus obtained 
from the data fit to the Landau free energy, Eq.Q. The full 
circles are for 64 Ni, the full triangles are for 70 Zn and the 
squares for 64 Zn projectiles impinging on various targets, see 
text. 



and c as a function of the compound nucleus m cn . As 
we see there is some dependence which may reflect dif- 
ferences in temperature. However, we note that the error 
bars and fluctuations are large which may also indicate 
important secondary decay effects. Thus is not too so 
easy to draw definite conclusions. 
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FIG. 7: The parameters a,b and c versus (I/A) of the com- 
pound nucleus obtained from the data fit to the Landau free 
energy, Eq.©. The symbols are like in figure [6] 

Given the information on the parameters of the 
Landau free energy contained in figures |5] and [7] we can 



discuss some features regarding the NEOS. In particular 
for each reaction system we can estimate F/T when 
H/T = 0. In figure [8] we plot this quantity versus m 
of the fragments for various reactions. The curves do 
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FIG. 8: F/T (H/T=0) versus m of the fragments obtained 
from the a, b and c parameters fit to the Landau free energy, 
Eq.(3). Results for all experimentally investigated reactions 
are displayed. 

not differ much suggesting that temperatures are quite 
similar. The fits exhibit curvature near m = ±0.4 
which may suggest the presence of additional minima 
at larger absolute values of m. This could indicate 
either a first order phase transition or superheating 
(see below). The lack of data at very large m makes 
it difficult to constrain the fit. However, we can study 
other situations of particular physical interest which 
arise when the relationships among the parameters a, b 
and c are constrained. |2H4|: 

We have considered four such cases as follows: 
1) superheating. This case corresponds to b = —\J\ac 
and gives two minima at m ^ and is plotted in fig- 
ure [5] for the 64 Ni+ 232 Th system with a short dashed 
line. These are not absolute minima, which occur only at 
m — 0, and they correspond to metastable states. They 
might be observed in high quality data for collisions of 
more neutron rich or proton rich systems making a hot 
source with m source rs ±0.4. In fact if the system could 
be gently brought to the right temperature T s , with the 
correct isotopic composition, it might stay in the mini- 
mum, i.e. more fragments of that m should appear; 
2) line of first order phase transition. This corresponds 
to the condition b = —y/16ac/3 at a temperature T3, 
which, if imposed on the fit of the free energy, results in 
the dashed-dotted-dotted-dotted line of FigfS] This fit is 
of similar quality to the previous cases. Now the minima 
are atm« 0.6, i.e. for more neutron rich fragments due 
to the fact that H/T 7^ 0. This suggests that in this 
situation we might produce a large number of neutron 
rich fragments. However, most of those fragments are 
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probably unstable, thus coincidence measurements may 
be required to determine their yields. Of course this fea- 
ture should become important in neutron rich stars. 

3) first order phase transition. Corresponds to the case 
a = and determines the 'critical temperature' Tq where 
the minimum at to — disappears and only the ones at 
to =/= survive. This case is excluded by our present data. 
However, the fit in Fig. ([5]) suggests an intermediate sit- 
uation between this and case (2) above. 

4) line of second order phase transition, tri-critical point. 
Corresponds to a = and b > (T = T c ). When b = as 
well we have a tricritical point (T = T$ c ), i.e. the point 
where the line of first order phase transition terminates 
into a second order phase transition. This case is also 
excluded by our data. 

We can extrapolate the cases discussed above to H/T = 
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FIG. 9: F/T (H/T=0) versus m of the fragments obtained 
from the a, b and c parameter fit to the Landau free energy, 
Eq.© for 64 Ni+ 232 Th. The four curves correspond to va- 
por(full line), superheating(short dashed), line of first order 
transition(3 critical dashed-dotted-dotted-dotted line), first 
order phase transitioning dashed line), see text. 

as was done for FigEl In Fig|3]we plot F/T (H/T=0) 
(extrapolated from the data) vs. to. Purists will not call 
this the EOS but reserve that for the pressure vs. to case 
(that we discuss below). Since H/T is zero, the curves 
are symmetric with respect to to. We see in the plot: 
vapor (dashed-dotted-dotted line) T > T s , superheating 
(T = T s ) (dotted line), a point in the line of a first or- 
der phase transition (T = T 3 ) (dashed-dotted line) which 
displays three equal minima at to and m±, see EqJSl 
the experimental data fit line (full line)T < T3. We have 
also added the case a — which should be obtained at 
T = To, where the minimum at to = becomes a maxi- 
mum. There is a series of cases not displayed in the fig- 
ure, corresponding to the temperatures between To and 
T3 where to = is still a minimum but not an absolute 
minimum. This corresponds to supercooling and might 
be observed in gentle collisions of ./V = Z nuclei similarly 
to the superheating case. 

The features in Fig[S]are reminiscent of the superfluid A 



transition observed as some 3 He is added to 4 He 0] . Pure 
4 He has a critical temperature of 2.18 K. The critical 
temperature for the second-order transitions decreases 
with increasing 3 He concentration until at temperature, 
T = 0.867K, a first-order transition appears. This point 
is known as the tri-critical point for this system. In a 
similar fashion, a nucleus, which can undergo a liquid- 
gas phase transition, should be influenced by the different 
neutron to proton concentrations. Thus the discontinuity 
observed in Fig[5] (to = 0) could be a signature for a tri- 
critical point as in the 4 He- 3 He case. We believe that 
our data, analyzed in terms of the the Landau 0(to 6 ) 
free energy, suggest such a feature but are not sufficient 
to clearly demonstrate this. Some other work [2^, [23 |. 
also suggests that a line of critical points might be found 
away from its 'canonical' position, i.e. at the end of a 
first-order phase transition and, for small systems, even 
extending into the coexistence region. 



III. CRITICAL EXPONENTS 

In the fits discussed above the parameters a, b and c 
were left free since we do not have any particular values 
to fix the scale. Nevertheless, we saw in figure [5] that the 
free energy (H/T=0) looks very similar for the different 
systems. Thus the values of the fitting parameters are 
similar apart from a scaling factor. We can avoid unnec- 
essary factors by defining suitable dimensionless quanti- 
ties. This can be accomplished by looking at the solutions 
of the minima of the free energy, cf. Eqf5] In particular 
from the value at the minimum, to + we can define the 
following quantities (b 7^ 0): 



4ac 



(7) 



Recalling that a is related to the distance from the crit- 
ical temperature while b and c should only depend on 
density [2J, we deduce that x is a measure of the dis- 
tance T — T3 from the critical temperature in a suitable 
dimensionless fashion. Similarly we can define a reduced 
order parameter from EqJS] 



2cm 2 



" |6| 

Thus Eqf5]can be rewritten as: 



y = l + Vl^ 



(8) 



(9) 



Near the critical point we know that the order parameter 
has a singular part that behaves in a power law fashion, 
thus we can define the singular part as: 

M = ±^l = ±{l-xf'\ (10) 

defining the temperature 'distance' from the critical 
point, |t| = |1 — x\, immediately gives the value of a 
critical exponent: j3 = |. This exponent is very close to 
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FIG. 10: Order parameter versus reduced temperature for 
all studied systems. The dashed line is given by Eq. (|10[) . the 
vertical line indicates the critical temperature T3. To the right 
of this line the system is in a superheated state. Supercooling 
occurs on the left of the vertical line and M=0 Q. 



the accepted experimental value as is well known in the 
0(m 6 ) Landau theory Q. In Fig. fTU)) the experimental 
values of M and x obtained within the Landau theory are 
plotted together with the equilibrium condition given by 
Eq.(fTU|). Supercooling and Superheating regions, as dis- 
cussed in the previous section, can be identified as well 

!• 

As is the case for macroscopic systems we can now 
'turn ' the external field H on and off. In our case this is 
done with a suitable choice of the colliding systems. In 
this way we can study the 'EOS' at the critical point by 
turning on H: 



M = H 1/s , 



(11) 



which defines the critical exponent 5. In the Landau 
theory this exponent can be determined at the critical 
point where a — b = 0. From equation Q we easily get 
<5 = 5, which is the accepted value for such a critical ex- 
ponent 0. In order to exactly determine this exponent 
we need to bring the system to the critical point. This 
does not appear to be the case for our data as we saw 
in figure [TU] Nevertheless a plot of the order parameter 
versus H should display a power law behavior as it is well 
known in macroscopic systems @. A precise determina- 
tion of the critical exponent requires the knowledge of the 
temperature T both above and below the critical point. 
This is feasible but requires precise experimental data. 
^From Eq. © assuming the only minimum is at m=0 we 
get 



V 



H/T 



(12) 



theless we can study the behavior close to the critical 
point by suitably defining scaling forms 0: = -S^ 

vs. Tjpr- These quantities are plotted in figure ITT1 and 
compared to magnetization data for nickel metal. The 
scaled magnetization is plotted versus the scaled exter- 
nal magnetic field [2[ . The nuclear data have been shifted 
in the region near the crossing of data above and below 
the critical temperature where we expect our data to be, 
see FiglTni Of course it is not possible at this stage to di- 
rectly compare to the macroscopic data since we have no 
information for the absolute values of the temperatures. 
Furthermore the role of the density (or pressure) is not 
clear since we expect that the parameter a (or equiva- 
lently x) depends on the 'distance' from the critical tem- 
perature and critical pressure. These quantities could 
however be obtained in 4-7T experiments where charges, 
masses and their velocities are carefully determined. 




The temperature (a) dependence of the order parameter 
shows that we are away from the critical point. Never- 



FIG. 11: Scaling form for magnetization M vs. external field 
for nickel |2|, open symbols. The corresponding quantities 
for nuclei normalized to the metal case are given by the full 
symbols. 

Once we have derived the 'reduced' parameters of the 
Landau 0(m 6 ) theory, we can write a 'reduced' free en- 
ergy as (b ^ 0): 

/ 1 2 ||4|2 2 6 ^ /1 n\ 

— = -xm - \z\m + -z to - —to, (13) 

where: = |§|J, Z = t§t, ^ = fifr. These quanti- 
ties together with the temperature, Eq.Q, and the re- 
duced order parameter y,of Eq.®, constitute the Landau 
0(m 6 ) theory in dimensionless form. It is instructive to 
study how these quantities change with the reaction sys- 
tem as we did in figures © and ©. In Fig. (IT2l we plot 
these normalized quantities vs. difference in neutron pro- 
ton concentration of the compound nucleus. Compare to 
figure [7J A feature worth noticing is the following, while 
the parameter a is decreasing with increasing (I/A) of the 
compound nucleus, the opposite holds for the parameter 



9 



0.8: 
0.6: 
0.4; 
0.2: 

1.5 



N 




*b 0.05 0.1 0.15 0.2 0.25 

(l/A) 

x 'cn 

FIG. 12: The parameters x, z and h/T versus (I/A) of the 
compound nucleus obtained from the data fit to the Landau 
free energy, Eq.(|3]). The symbols are like in figure [6] 



x which gives the 'distance' from the critical temperature, 
see Fig. (|12[) . This is very important since only normal- 
ized quantities should be used when inferring the prop- 
erties of the EOS (i.e. temperature, density etc.) near 
the critical point. 



IV. SYMMETRY AND PAIRING COMPARED 
TO THE COULOMB ENERGY. 

In the previous sections we have seen that the Coulomb 
energy might become important especially for large val- 
ues of the charges. We can now try to derive some qual- 
itative understanding of when and why Coulomb correc- 
tions might become important and might even hinder a 
possible phase transition. From the mass formula we can 
write the Coulomb energy for large Z as Q]: 

^£ = .77^U 2 / 3 = — {\-mfA 2 l\ (14) 
A A 2 4 v ; v ' 

which explicitly introduces the order parameter to in the 
Coulomb energy. We can define an 'effective' symmetry 
energy (per particle) as: 



E °" = (a sym + ^A 2 / 3 )m 2 ^ A^m + "-^A 2 ' 3 , 



A 



0.77 
~2~ 



0.77 

T 



(15) 

where the symmetry energy coefficient a sym = 25MeV. 
Ignoring for a moment density corrections we see that 
0(m 2 ) term should be affected by Coulomb corrections 



for large fragment mass numbers. Furthermore, a linear 
term in m is introduced which will then modify the 'ex- 
ternal' field even in collisions where the source m s = as 
we discussed in Fig.©. Finally there is a term not de- 
pendent on m that will destroy the scaling for large mass 
(charge) numbers. We should also notice that assuming a 
spherical expansion, at low densities the Coulomb energy 
will decrease as p 1 ^ 3 while contributions to the symmetry 
energy should depend both on p 2 / 3 reflecting the Fermi 
energy of the nuclei and on p, the latter coming from 
different n-p interactions. At low densities we would ex- 
pect Coulomb to be stronger than it appears to be in the 
data. This may be indicative that the fragments must be 
highly deformed, reducing the Coulomb energy. Coulomb 
corrections should become more important when to = 
for the detected fragment. We have plotted the yields of 
to = nuclei in Fig. ([2]) and pointed out that pairing ap- 
pears to be playing a role. From Eq.(fT5]) above we should 
expect that, if Coulomb is dominant for such fragments, 
the free energy should depend on A 2 ' 3 . In Fig. fllUl) (top 
panel) we plot F/T versus A for m = fragments. The 
expected dependence with mass number in the free en- 
ergy suggested from 'effective' symmetry energy Eq. (|15[) 
is not seen in the figure. Rather, a staggering between 
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FIG. 13: Free energy versus mass for m = isotopes for the 
70 Zn-|- 124 Sn system (top panel), the dashed line is a fit using 
Coulomb and pairing contributions. Free energy times S (see 
text-bottom panel) versus mass for m — isotopes. The lines 
are separate fits suggested by the Coulomb (dashed line) and 
pairing(full line) energy mass number dependence. 



odd-odd and even-even nuclei is clearly visible. 

To better clarify these arguments we can write the pair- 
ing energy from the mass formula as |l| : 



E P _ U _S 



(16) 



where 12MeV is the ground state pairing energy coeffi 
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cient and <5: 

(-1) N + (-l) z 
6= [ ' + } ■ (17) 

The suggested mass dependence from pairing, Eq |16[ is 
completely different from the Coulomb one when m = 0, 
see Eg mi 

Notice that it is the 5 factor from pairing that changes 
the sign of the contribution for odd-odd to even-even nu- 
clei. A combined fit to the data using Coulomb plus pair- 
ing contributions results in the dashed line in Fig.(fT3"|). 
top panel. The agreement with data is very good. If 
we multiply the pairing energy by the factor S we should 
get no discontinuities when plotting this quantity versus 
mass number. Similarly if the properties of the free en- 
ergy depend on the pairing term, as for the ground state 
case, then it should be a monotonic function of A after 
multiplying it by 5. In Fig. fTB")) (bottom panel), we plot 
the quantity £<5 versus mass number for the same sys- 
tem of Fig. (IT3l) (top panel). The fit using the pairing 
mass dependence is also good. The Coulomb mass de- 
pendence fails especially for small mass number. From 
the values of the fit, using the ground state coefficients 
we can derive a temperature for the Coulomb case of 
T = 9.2(-j^) 1 / 3 MeV where we have explicitly indicated 
a possible density correction. For the pairing case we get 
T = 6A5MeV. Notice that in this case we have not sug- 
gested any density correction since the fate of the pair- 
ing energy at low density and finite temperature is 'terra 
incognita'. When making a combined fit using pairing 
and Coulomb energy we get a good reproduction of the 
data (dashed line in Fig. (fT3j) .top panel). While the fit- 
ting value for pairing results in a 'temperature' T=5.13 
MeV, we get an increase of the Coulomb contribution to 
T=12.1 MeV. Assuming that pairing is independent on 
density, we could derive a density from the Coulomb re- 
sult. A simple calculation give = (6.45/9) 3 = 0.34 
which could be a reasonable indication of the density of 
the system when it breaks into fragments. 

In summary in this section we have shown that the role 
of the Coulomb energy appears to be rather reduced in 
the reactions analyzed in this paper. We expect it to be- 
come more important for large nuclei. On the other hand 
large nuclei have smaller symmetry and pairing energies 
per nucleon, thus a precise determination of the EOS can 
be obtained from measurements of isotopes having rela- 
tively small masses. 



V. DYNAMICS OF THE PHASE TRANSITION 

As we have seen we have been able to discuss some ob- 
servables in the fragmentation of nuclei using a language 
common to macroscopic systems undergoing a phase 
transition. In the nuclear case we have a finite system 
composed at most of hundreds of particles which evolves 
in time under the influence of a long range Coulomb force. 



This poses many questions on why techniques of statisti- 
cal mechanics should apply in such evolving nuclear sys- 
tems. This also offers the possibility of dealing with sta- 
tistical mechanics of open systems and the problem of 
extending the description of a phase transition to such a 
system. 

We start by observing that even though we are deal- 
ing with a dynamical system, the order parameter de- 
fined in this work, m, is confined between -1 and +1. In 
this sense we have a somewhat 'closed' system. Also the 
density at which the transition occurs should be smaller 
then normal density and thus Coulomb effects are re- 
duced. However, if we deal with larger sources, such as 
in U+U collisions the phase transition might be washed 
out by the strong Coulomb field. We expect our current 
considerations to be valid for small sources only. 

From statistical mechanics we know that in a first order 
phase transition |2j a small seed increases in size depend- 
ing upon the surface tension at a given T and density p. 
If the pressure of the surrounding matter is smaller than 
the internal pressure of the drop, the drop will grow by 
capturing surrounding matter. On the other hand if the 
opposite is true then the drop will decrease in size to bal- 
ance the external pressure. The entire process is driven 
by surface tension. Drops of a given size will survive 
only when their internal pressure balances the external 
pressure. If the system is at a very low density the in- 
teraction between different parts might take a relatively 
long time. In these conditions a big nuclear drop whose 
internal pressure is larger than that of the surroundings 
could be considered to be a nucleus which is evaporat- 
ing particles in order to balance the external (zero in the 
case of an isolated nucleus) pressure. If we accept this 
picture than the evaporation step is part of the dynamics 
of the phase transition. Thus a very low density system 
might be thought of as many isolated drops evaporat- 
ing particles and reaching their equilibrium conditions 
before they collide with other parts of the system or as 
small fragments being evaporated by other drops. In a 
finite system this does not happen, but we might think of 
a process where at some point the finite system becomes 
unconfined and an infinite system is approximated by 
an infinite number of repetitions or 'events'. Of course 
in a statistically equilibrated system we know that time 
averages and event averages are the same. Here we are 
extending this concept to finite systems where only event 
averages can be used. A major question here is whether 
the properties of the phase transition are decided very 
early, i.e. when the system 'enters' the instability re- 
gion. As we said above if we have an infinite system at 
a very low density undergoing a first order phase tran- 
sition, then the drops can explode, evaporate, and fuse 
with other particles over a very long time. Our finite 
system might behave similarly but without the fusion at 
later times. If this were the case than the detected frag- 
ments carry all the information of the phase transition, 
if not then we need to reconstruct the primary fragment 
distributions coming out of the instability region. 



11 



< 

o" 



2 

-2 





\ / 


• 

• • 


\ ? 


-* ' ' ' r 

\ . • f 


A . . f 



-1 -0.5 0.5 1 -1 -0.5 0.5 1 
(N-Z)/A (N-Z)/A 



FIG. 14: Free energy vs. time in AMD calculations (see text) 
for 64 Zn+ 112 Sn system at 40MeV/A and central collisions, i.e. 
impact parameter less than 3 fm. Different picture correspond 
to T=200, 300, 500, 1000, 1500 and 2000 fm/c respectively. 



We can try to clarify some of these questions by means 
of microscopic models such as Antisymmetrized Molec- 
ular Dynamics (AMD) or similar approaches where the 
time evolution of the system is followed [22| . However, we 
have to stress that in such microscopic models some as- 
sumptions are made in order to recognize the fragments 
at particular times during the time evolution. In sim- 
pler approaches, fragments are recognized if particles are 
close in coordinate space (of the order of the range of the 
attractive nuclear forces) Q. In such a case the recog- 
nized fragments are 'excited' and they evolve in time until 
a final state is reached after a long time of the order of 
thousands of fm/c. A more refined approach for fragment 
recognitions is given by defining clusters when its com- 
ponents are within a given distance in phase space. The 
naive expectation would be that in this case we should 
recognize fragments earlier than the previous case and 
this is the method that we will adopt here for simplic- 
ity following ref . [22[ . In an ambitious approach [l(| the 
claim is that fragments are recognized very early during 
the time evolution, of the order of tens of fm/c, if one 
searches for particles connected in phase space to form 
fragments and minimize the energy. This case probably 
corresponds to minimizing the entropy of the fragment- 
ing system. If this last picture will hold true then a 
picture of an infinite system at low density will be equiv- 
alent to an 'infinite' repetition of events. Finally in all 
the considerations above we have to add the necessary 
and interesting complication that we have a mixture and 
not a single fluid, thus we can have more situations to 
explore than discussed in the previous sections and we 
can 'turn on and off' an external field as well. 

We have performed AMD calculations for the same 
systems investigated experimentally. After some time t, 
fragments are separated enough in phase space that we 
can recognize within a simple phase space coalescence ap- 
proach as discussed in [22|. In this way we can define a 
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FIG. 15: Fit parameters a, b and c vs. time (see text) for the 
same system of Fig ll4l Solid circles refer to AMD calculations 
while the open symbol in the M vs. x plot is the experimental 
value for this system. 



yield at a given time and from this derive the free energy 
exactly as we did with the experimental data. Charac- 
teristic results for the free energy versus time is given in 
Fig. ([l~4")) together with a Landau 0(m 6 ) fit. Some time 
evolution is observed. Using a more sophisticated frag- 
ment recognition approach [10] might even decrease the 
time over which this evolution occurs. We can study the 
time evolution in more detail by plotting the variables 
a, b and c. M defined in the previous sections versus 
time. The results of the fits to the free energy at dif- 
ferent times is given in Fig. (fT5|) . While the quantities 
a, b and H/T change somewhat during the time evolu- 
tion, smaller changes are observed in the time evolution 
of normalized quantities, x, z and h/T. Nevertheless the 
time evolution of the fitting parameters influences the 
time evolution of the order parameter M versus reduced 
temperature x as seen in the bottom right of Fig.tfPoj). 
It is very interesting to see that in these units the sys- 
tem is initially very hot (superheated) and cools down 
when coming to equilibrium below the critical tempera- 
ture. The final result is very close to the observed values 
given by the open points. 

Thus in this model most qualitative features of the 
phase transition are decided very early during the time 
evolution. This might correspond to an entropy satura- 
tion early during the evolution. However, different mod- 
els and fragment recognition approaches might change 
the picture somewhat. 
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VI. EQUATION OF STATE 

Once we know the free energy (at least in some cases) 
we can calculate the NEOS by means of the Fisher 
model [2l| . Since we do not have at present experimental 
information on the density p, temperature T and pres- 
sure P of the system we can only estimate the 'reduced 
pressure' pH ]: 



(18) 



where Mj are moments of the mass distribution given by: 



Mi, 



A 



A k Y{A,m) =Y J2 AkA ~ 



-F/T(m)A. 



0,...n. (19) 



Notice that the quantities above are now dependent on 
the order parameter m. From the knowledge of F/T 
(H/T = 0) from the previous section we can easily cal- 
culate the reduced pressure near the critical point. In 
particular given the simple expression for the moments 
we can also derive some analytical formulas following [25[ : 



P 



(m) 



3.072|F/T| 4 / 3 + 1.417 - 3.631|F/T| -t 
-4.086IF/TIV3 + 3.631 + 0.966|F/T| 



which gives at the critical point a critical compressibility 
factor {F/T = 0): = ££g =0.39. 
This value is essentially that derived from the Van-der- 
Waals gas equation but is well above the values observed 
for real gases. Using the relations above we can calculate 
the NEOS for the situations illustrated in FigEl The 
results are displayed in FigOJj]where the reduced pressure 
is plotted versus m for vaporization, superheating and 
first order phase transitions on the tri-critical line. Notice 
that there is not a large difference between the first two 
cases, while the last case displays two critical points (a 
third one is on the negative m axis). 
We have seen in FigJ5]that N = Z nuclei display a power 
law. We can also estimate the critical reduced pressure 
for this case noticing that the sums in Eq. (fT9l) above must 
be restricted to A = 2Z nuclei. This leads to a critical 
compressibility factor 



0.20 which is a value closer 



to that estimated from other multi-fragmentation studies 
before [13 . 

We can compare our analytical result given in Eq. (|20[) 
with the numerical values obtained above. This is dis- 
played in Fig. (|17|) and we see that the numerical approx- 
imation is especially good near the critical point (s) as 
expected. If, from detailed comparison to experimental 
data, we are able to extract the temperature and pres- 
sure dependence of the parameters entering the Landau 
free energy, then Eq. (fT8)) would be the Nuclear equation 
of state near a critical point. From the actual data at 
our disposal we can only estimate the behavior of the re- 
duced pressure as function of the order parameter m. On 
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FIG. 16: Reduced pressure versus m of the fragments ob- 
tained from the a, b and c parameters fit to the Landau free 
energy, Eq.((3]) for the 64 Ni+ 232 Th. The curves correspond 
to vapor(open circles), superheating(open squares), first or- 
der (3 critical line-solid stars), see text. The solid circle is for 
N = Z nuclei at the critical point. 
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FIG. 17: Comparison to the analytical result, Eq. (|20[) (solid 
circles) for a first order phase transition. 



similar ground we can define a reduced compressibility as 
XP T(m) = ^. (21) 



Its behavior is displayed in Fig. (fl8l) for the cases out- 
lined above. Divergences near the critical point(s) are 
obtained. 



VII. CONCLUSIONS 

In conclusion, in this paper we have presented and 
discussed experimental evidence for the observation of 
a quantum phase transition in nuclei, driven by the neu- 
tron/proton asymmetry. Using the Landau approach, we 
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FIG. 18: Reduced compressibility versus m of the fragments 
as in Fig.|l6j) 

have derived the free energies for our systems and found 
that they are consistent with the existence of a line of 
first-order phase transitions terminating at a point where 
the system undergoes a second-order transition. The 
properties of the critical point depend on the symmetry. 
This is analogous to the well known superfluid A transi- 
tion in 3 He- 4 He mixtures. We suggest that a tricritical 
point, observed in 3 Hc- 4 He systems may also be observ- 
able in fragmenting nuclei. These features call for fur- 
ther vigorous experimental investigation using high per- 
formance detector systems with excellent isotopic iden- 
tification capabilities. Extension of these investigations 
to much larger asymmetries should be feasible as more 
exotic radioactive beams become available in the appro- 
priate energy range. 

It is important to stress that the observables dis- 
cussed here represents only necessary conditions for a 
critical behavior. A definite proof of a phase transition 
and a tricritical point could be given by a precise 
determination of yields of fragments whose m f=a ±0.5, 



i.e. very unstable nuclei which, most probably, decay 
before reaching the detectors. Thus fragment-particle 
correlation measurements for exotic primary fragments 
such as A Li, 5 Be (proton rich) or extremely neutron 
rich 10 He are needed. More generally, such correlation 
experiments can also shed light on the effects of sec- 
ondary decay on the fragment observables. This remains 
a key question in many equation of state studies and 
model calculations differ in their assessment of these 
effects d, 0] ■ Higher quality data over a wider range of 
beam energies and colliding systems should also help 
in clarifying the role of other energy terms, such as 
surface, Coulomb etc., which are important at lower 
excitation energies. In particular the role of pairing and 
the possibility of Bose-Einstein condensation, should be 
more deeply investigated. Our data for I — fragments 
already show that pairing is important. This might be 
due to its importance during the phase transition or to 
its role during secondary decay of the excited primary 
fragments. Exploration of quantum phase transitions in 
nuclei is important to our understanding of the nuclear 
equation of state and can have a significant impact in 
nuclear astrophysics, helping to clarify the evolution of 
massive stars, supernovae explosions and neutron star 
formation. 
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